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Abstract 

We consider a model which in a certain limit reduces to the large N N = 1 su- 
persymmetric SU(N) gauge theory without matter. The gaugino condensate 
in this model is controlled by the dynamics of an additional singlet superfield. 
Using this model we explicitly construct BPS domain walls arising due to the 
chiral symmetry breaking. In particular, in the large N limit we obtain the 
exact shapes of the domain walls corresponding to solitons, and also of the 
domain walls interpreted as D-branes on which the SQCD string can end, 
whose existence was previously argued by Witten in the context of the large 
N SQCD. We also discuss various points which appear to support the consis- 
tency of the D-brane interpretation for these domain walls within the SQCD 
string context. 
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I. INTRODUCTION 



One of the longstanding problems in high energy physics is to gain analytical control over 
strongly coupled gauge theories, and, in particular, QCD. Many fascinating aspects of QCD, 
such as confinement, mass gap and chiral symmetry breaking, still lack desired quantitative 
description, albeit over the years much qualitative (and numerical) insight has been gained 
into these phenomena. The closest supersymmetric relative of QCD, namely, pure TV = 1 
super symmetric SU(N) gauge theory a priori offers a possibility of simplifications (due to 
super symmetry), yet it shares most of the qualitative features of QCD. Thus, SQCD without 
matter is believed to be a confining theory with a mass gap and chiral symmetry breaking. 
It is therefore reasonable to suspect that a better understanding of pure SQCD might shed 
light on these phenomena in a more general context. 

It is not unreasonable to believe that confinement and chiral symmetry breaking might 
be intimately related JTJ. If so, understanding chiral symmetry breaking in SQCD might 
shed some light on the confining properties of (super)glue. Chiral symmetry breaking in 
SQCD is due to gaugino condensation. Thus, the composite operator AA acquires a non- 
zero expectation value in strongly coupled SQCD: (AA)fc = a(N) exp(2iiik/N)A 3 , where A is 
the dynamically generated scale of SQCD, and a(N) is an overall factor which depends on 
the subtraction scheme. There are A" inequivalent vacua labeled by k = 0, . . . ,N — 1. This 
implies that there should exist domain walls connecting pairs of distinct vacua. In |2| these 
domain walls were argued to be BPS saturated. The BPS property of the SQCD domain 
walls suggests that we might be able to gain some analytical control over these solutions, 
and, therefore, over chiral symmetry breaking. 

A seemingly unrelated subject is 't Hooft's large N limit ||. In it was noticed that 
in this limit the gauge theory diagrams organize themselves in terms of Riemann surfaces 
characterized by the numbers of boundaries and handles, and addition of each extra handle 
on the surface corresponds to suppression by a factor of 1/N 2 . This observation has led 
't Hooft to speculate that the large A" gauge theory might be described by some kind of 
string theory. Then the large A" expansion of gauge theories would be mapped to the string 
expansion in terms of the properly weighted world-sheets of various topologies. 

In the context of weakly coupled gauge theories (i.e., when the effective gauge coupling 
A = NgyM is fixed at a value A < 1) the first concrete realization of this idea was given in || 
in the context of three dimensional Chern-Simons gauge theory where the boundaries of the 
string world-sheet are "topological" D-branes. More recently, string expansion was shown to 
precisely reproduce 't Hooft's large A" expansion for certain four dimensional gauge theories 
||[]. Here the boundaries are Type IIB D3-branes, and the (open) QCD string is simply the 
open string stretched between these D-branes. 

As to the case of strongly coupled gauge theories, the story here appears to be much 
more complicated. In particular, it is unclear what is the string theory that governs the 
dynamics of the QCD string which is expected to arise as an effective description in the 
large A" limit. This string theory is expected to be non-critical M, which makes it difficult 
to study. Nonetheless, in || Witten has made a remarkable observation that there might 



This was subsequently generalized to include unoriented world-sheets in ||. 
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be a connection between SQCD string and domain walls in the large N limit. In particular, 
based on the assumption that the SQCD domain walls are BPS saturated, Witten argued 
that in the large N limit the domain walls connecting two vacua labeled by k and k! = k + 1 
appear to be objects that are not solitons from the SQCD string viewpoint but rather look 
like D-branes || on which the SQCD string should be able to end. Such D-brane-like domain 
walls have recently appeared in [plj in the context of the AdS/CFT correspondence ||11|| . 

It is desirable to verify the above assertions explicitly since they have important implica- 
tions for SQCD. One of the key observations is that there exist D-branes on which the SQCD 
string can end ||. In particular, the chromo-electric flux contained in the SQCD string that 
ends on the SQCD D-brane spreads along the D-brane, but away from the D-brane it is 
trapped in the thin flux tube corresponding to the SQCD string [j3Ul2|. In fact, a test elec- 
tric charge can only come off a D-brane (in the transverse direction) if it is connected to it 
by the SQCD string. Another way of putting this statement is that the D-branes should 
support excitations that behave as if they transform in the fundamental of SU(N) albeit 
there are no such fields in the original super- Yang-Mills Lagrangian ||. Moreover, naively 
one expects that the pure superglue SQCD contains only closed strings. However, if large 
iV domain walls indeed behave as D-branes, then we must also include open strings in the 
corresponding (non-critical) string theory to be able to describe these domain walls. This, 
in a sense, is analogous to what happens in Type II string theory: perturbatively we only 
expect closed strings; however, to capture non-perturbative dynamics we must also include 
open strings which end on the corresponding D-branes. 

The purpose of this paper is to give more evidence supporting the D-brane picture for 
the large N SQCD domain walls. In particular, using an indirect construction, we obtain 
solutions for the large N SQCD domain walls which are BPS saturated. Also, it is important 
to check that the D-brane picture for the domain walls connecting two vacua labeled by k 
and k' = k + j with \j | = 1 is consistent with the fact that the domain walls connecting 
such vacua with \j\ ~ N are expected to behave as solitons of the SQCD string. Moreover, 
a priori it is far from being obvious that the domain walls connecting two vacua labeled by 
k and k' = k + j with \j\ ~ N a , < a < 1, are consistent with the effective SQCD string 
description. In particular, what is the interpretation of such domain walls in the context of 
the SQCD string? 

It is, however, difficult to explicitly construct these domain walls (or even check that they 
are BPS saturated) directly in the context of SQCD: the order parameter AA is a composite 
operator, and we are dealing with a strong coupling regime. Nonetheless, in this paper we 
use an approximating model which in a certain limit reduces to the large iV SQCDQ. In fact, 
using this approximating model we will be able to solve for the exact shapes of domain walls 
for \j\ = 1 (these domain walls are interpreted as SQCD D-branes), and also for \j\ ~ N 
(these domain walls are interpreted as solitons of the SQCD string). More precisely, we will 
obtain the corresponding exact solutions in a certain coordinate system which depends on 
the Kahler potential. The latter we will assume to be non-singular (in the regime we are 
interested in here). Moreover, it will become clear that these solutions are BPS saturated. 



2 More precisely, it appears to be in the same universality class as the large N SQCD meaning 
that the two theories are the same at low enough energies. 
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Also, we will argue that for all other values of \j\ the corresponding domain walls in the 
large N limit are described by \j\ D-branes separated by large distances. This provides an 
important check for the consistency of the D-brane interpretation and fits in the SQCD 
string picture. 

The exact results we obtain in this paper are partly due to the fact that the domain wall 
problem reduces to a two dimensional problem, namely, of finding soliton configurations in 
two dimensional M = 2 massive quantum field theories. These theories have been studied 
extensively (in fact, one of the original motivations to study them was precisely their rele- 
vance to domain walls in four dimensional M = 1 supersymmetric (gauge) theories), and the 
vast knowledge of their properties should not be overlooked when studying domain walls. 

The rest of this paper is organized as follows. In section II we review some well known 
facts about M = 1 supersymmetric QCD. In section III we give a general discussion of 
domain walls in M = 1 supersymmetric theories. In section IV we construct the approxi- 
mating model, and solve for the domain wall shapes in the large N SQCD. In particular, 
we explicitly construct the BPS solutions corresponding to the D-branes and solitons. We 
also identify the limit in which the approximating model reduces to the large N SQCD. In 
section V we discuss the zero modes localized on the domain walls. In particular, we identify 
the zero modes localized on the D-branes (which are described by a three dimensional field 
theory) with the U(l) gauge supermultiplet in three dimensions as required by the D-brane 
interpretation. In section VI we discuss various issues regarding the large N SQCD string 
and domain walls. 

II. PRELIMINARIES 

In this section we review some well known facts about M = 1 supersymmetric QCD. 
In particular, we briefly discuss gaugino condensation in SQCD without matter. We then 
consider domain walls in SQCD. In the subsequent sections we will argue that these domain 
walls are BPS saturated. 

A. Gaugino Condensation in SQCD 

Consider M = 1 supersymmetric QCD with SU(N) gauge group and no matter. Let 
A be the dynamically generated scale. More precisely, if g(/i) <C 1 is the Yang-Mill gauge 
coupling constant at some energy scale /i, then A = /iexp(— 87i 2 /b g 2 ([i)) is the scale at 
which the one-loop renormalized gauge coupling blows up. (Here b = 3N is the one-loop 
/3-function coefficient.) 

In this theory there is a non-perturbative superpotential due to gaugino condensation: 

W = N(XX) . (1) 

The gaugino condensate is given by 

(AA) = a(N) exp(2mk/N)A 3 . (2) 

Here k = 0, . . . , N — 1, that is, there are N inequivalent vacua corresponding to N different 
phases for the gaugino condensate. This is due to spontaneous breaking of Zjy chiral symme- 
try. The factor a(N) (which a priori can depend on N) is a positive number, and it depends 
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on the subtraction scheme. However, the N dependence of a(N) becomes important in the 
large N limit, in particular, consistency of the SQCD string interpretation for the large N 
domain walls suggests that a(N) ~ iV [0]. This also naturally follows from the fact that 
the gaugino condensate (AA) involves a trace over the gauge indices, and the resulting vev 
should contain the corresponding second Casimir factor which in this case is simply N. 

The superpotential in ([]]) can be derived as follows ||14j| . Consider M = 1 supersymmetric 
gauge theory with SU(N) gauge group and N flavors Qi, Qi = 1,...,N). The gauge 
invariant degrees of freedom are mesons = QiQ\ and baryons B = e tl '" tN Qi 1 ■ ■ -Qi N 
and B = ej 1 ..rj N Q^ x • • ■ Q^ N ■ The classical moduli space in this theory receives quantum 
corrections which can be accounted for via the following superpotential flH|p|: 



W = A (det(M) -BB- A 2N ) . (3) 

Here A is the dynamically generated scale of the theory, and A is a Lagrange multiplier 
(AA 2N = W a W a is the "glue-ball" field). To obtain SQCD without matter, let us add a 
tree-level superpotential of the following form: 

W tree = mTr(M) . (4) 

This superpotential corresponds to giving mass m to all iV flavors. For m > A, at energy 
scales below m we can integrate the heavy flavors out which leaves us with SQCD without 
matter. The superpotential now reads 

W = Nexp(2nik/N)mA 2 . (5) 

The scale matching implies that 

mA 2 = a(N)A 3 . (6) 

That is, the scale matching depends on the threshold corrections (due to the heavy quarks 
whose mass is of order m), and, therefore, on the subtraction scheme. Thus, we arrive at 
(§■ 



B. Domain Walls in SQCD 

According to (jl]) and (@), there are inequivalent vacua in SU(N) SQCD without 
matter. We therefore expect that there exist domain walls separating these different vacua 
00. On dimensional grounds we expect the tension of such a wall to be of order A 3 . It 
is, however, possible to compute the tension exactly in terms of the gaugino condensate 



3 Here we note that this superpotential itself depends on the subtraction scheme. In particular, 
A is identified with the one-loop scale only up to a multiplicative constant that depends on the 
subtraction scheme. This constant equals 1 in the DR scheme. 

4 For other related works, see, e.g., [15]. 
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provided that they are BPS saturated 00. Consider the domain wall separating the vacua 
with (XX) k = a(N)exp(27rik/N)A 3 and (XX) v = a(N) exp(27rik' /N)A 3 . The tension of the 
wall (provided that it is BPS saturated) is nothing but the value of the central charge Qkk' 
(in the corresponding central extension of the M = 1 superalgebra) . The central charge is 
proportional to (the absolute value of) the difference between the values of the superpotential 
in the two vacua: 



It is, however, difficult to explicitly construct these domain walls (or even check that they 
are BPS saturated) directly in the context of SQCD: the order parameter AA is a composite 
operator, and we are dealing with a strong coupling regime. Nonetheless, in the subsequent 
sections we will see that using an indirect construction (i.e., by using a certain trick) it is 
possible to gain some insight into the structure of these domain walls. In particular, we will 
argue that these domain walls are indeed BPS saturated. Moreover, in the large N limit we 
will even be able to obtain the exact shapes of these domain walls. (More precisely, we will 
do so in a certain coordinate system which depends on the Kahler potential.) 

The presence of BPS saturated domain walls in SQCD has important implications ||. 
Strongly coupled SQCD in the large N limit is believed to be described by a string theory 
with the string coupling A s ~ 1/N ||. Extended solitons in this string theory are expected 
to have tension which goes as 1/A 2 ~ iV 2 . It is not difficult to see that not all of the 
above domain walls appear to be solitons from the SQCD string viewpoint. Thus, consider 
"elementary" domain walls (with k' = k+1). In the large N limit their tension goes as a(N). 
For these objects to be SQCD solitons we therefore would have to assume that a(N) ~ N 2 . 
Let us now consider the walls with \k' — k\ ~ N. The wall tension in these cases goes as 
Na(N) ~ iV 3 ~ 1/A 3 (provided that a(N) ~ N 2 ). This result is (at least) strange since we 
do not expect such dependence on A s for any BPS objects in the SQCD string context. 

To avoid the above problem with the \k' — k\ ~ N domain walls, we can assume that 
a(N) ~ N. In this case these objects have tension that goes as N 2 , and they can be 
interpreted as solitons of the SQCD string. Now, however, we see that the "elementary" 
domain walls have tension which goes as iV ~ 1/A S . These cannot therefore be interpreted 
as solitons. In it was suggested that such domain walls can be viewed as D-branes [[| in 
the SQCD string context. Open SQCD strings then can end on these D-branes ||. 

Here we would like to point out another important property of the domain walls in 
SQCD. Thus, consider the domain walls with \k' — k\ ~ N a , where < a < 1. If we 
assume that a(N) ~ N, then for these domain walls we naively get the tension that goes 
as N 1+a ~ l/Xl +a . Such states would be difficult to interpret in the context of the large N 
SQCD string. In the following we will argue that in the large iV limit these domain walls 



5 This was argued to be the case for SU(2) in Q. 
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(7) 



m - w k > 



We therefore have the following tension for the domain wall: 




(8) 
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actually correspond to \k' — k\ "elementary" domain walls separated by large distances. This, 
in particular, avoids the problem with interpreting these domain walls in the context of the 
large N SQCD string. 

Before we conclude this section, we note that if we chose a(N) to be independent of N 
we would run into a puzzle ||: the \k' — k\ ~ N domain walls would look like D-branes, but 
the "elementary" domain walls would have tension independent of the string coupling (and, 
therefore, would have to be considered as fundamental, that is, on the equal footing with 
the SQCD strings). 



III. DOMAIN WALLS IN N = 1 THEORIES 

In this section we provide some details concerning generic properties of domain walls in 
M = 1 supersymmetric theories. In particular, we discuss the relation between the domain 
walls in four dimensional M = 1 theories and solitons in (massive) M = 2 quantum field 
theories in two dimensions. 

Consider an TV = 1 supersymmetric theory with onef] chiral superfield X. Let W(X) be 
the superpotential in this theory such that the F-flatness condition Wx = has a discrete 
setQ of solutions X = X a , a = 1, . . . , N. 

In such a theory we expect presence of domain walls separating inequivalent vacua X = 
X a and X = Xf,, a ^ b. More precisely, let z be the spatial coordinate transverse to the 
wall. Then asymptotically we have X(z) —>■ X a as z — > — oo, and X(z) — > Xf, as z —>■ +oo. 
Such domain walls may or may not be BPS saturated. 

Note that the problem at hand is really a two dimensional problem as the coordinates 
x, y along the wall play no role in the discussion. Thus, our problem effectively reduces to 
that of finding BPS solitons in two dimensional massive M = 2 quantum field theories. More 
concretely, we are led to finding configurations of minimum energy (interpolating between 
a pair of distinct vacua), which are solitonic solutions of the Landau-Ginsburg equations. 



For a general Landau-Ginsburg theory |16] the energy of a time-independent configuration 
is given by 



l r + °° - / - |2 



oo 



E ab = - dz{g\d z X\ 2 + V) , (9) 

Z J—c 



where V = g~ 1 \Wx\ 2 is the scalar potential of the theory. Here g(X,X*) is the Kahler 
metric which is related to the Kahler potential /C via g = ICxx*- Note that E ab can be 
rewritten as follows: 

E ab = - dzg- 1 \gd z X* - exp(2 7 )Wx| 2 + Re (exp(z 7 )(W f) - W Q )) , (10) 



6 Generalization of the following discussion to cases with multiple superfields is completely 
straightforward. 

7 For the sake of simplicity we will assume that these vacua are non-degenerate. 
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where exp(«7) is an arbitrary constant phase. Note that E ab is independent of 7. For the 
choice 7 = ■yb a we obtain the bound E ab > \W b — W a \ which is the BPS bound. Here 

expH7fea) = w^m • ( } 

The BPS solutions (for which E ab = |W& — W a |) are those that satisfy the following equation: 

gd z X* = exp(t lba )W x (12) 
subject to the boundary conditions X(z) — > X a)b as z — ► =Foo. Going back to the domain 



walls in four dimensions, we have the exact same BPS equation (12), and the tension of the 
domain wall is given by 

T ab = -^\W a -W b \ . (13) 

Typically, one does not know the exact form of the Kahler metric g. However, for the 
case of a single superfield X we are considering here one can still make certain statements 
about the corresponding soliton solutions. Thus, let X be a function of a new coordinate z' 
such that it satisfies the following equation 

d z ,X* = exp(* 7ba )>Vx (14) 

with the boundary conditions X(z') — > X a<b as z' — > =1=00. Suppose we are able to find the 
corresponding solution for X(z'). Next, consider the following change of variables: 

d z! z(z')=g(X(z'),X*(z')) , (15) 

where X(z') is the corresponding solution. Note that if we express X as a function of z, 
then X(z) will satisfy the original BPS equation (|P2"|) with the corresponding boundary 
conditions. The change of variables ( |l~5f) is simply a diffeomorphism which is one-to-one as 
long as the Kahler metric g is non-singular. Throughout this paper (when discussing domain 
walls in four dimensional M = 1 theories or the corresponding solitons in two dimensional 
TV = 2 theories) we will assume that the corresponding Kahler metric g is indeed non- 
singular. Moreover, we will always work in the coordinate system parametrized by z', but 
we will drop the prime for the sake of simplicity of the corresponding expressions. (This is 
effectively equivalent to the case where g — 1.) Note that as long as the Kahler metric g is 
non-singular, the solution of fllil) (with the appropriate boundary conditions) implies that 
the corresponding solution of (O) exists. Moreover, this solution is BPS saturated. 



The solitons in Landau- Ginsburg theories have been studied in detail [^-|T8[. Many 



of these models are integrable. The general conditions for existence of BPS solutions were 



formulated in [18]. In this paper we will mainly be interested (for the reasons that will 
become clear in the next section) in the so called models for which the superpotential 
is given by 

X N+1 

W = A--_ . (16) 
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This model was studied in detail in 
model with X = = exp(2irik/N), k 
the vacua Xj, and X^ has the energy 



In particular, there are N distinct vacua in this 
0, . . . , N — 1. The soliton interpolating between 



E, 



kk' 



2N 
N + l 



sm 



'ir(k - k') 
N 



(17) 



All of these solitons are BPS saturated. (Here, as we have already pointed out, we take 
(7 = 1.) Here we also mention that the (purely elastic) S- matrix is exactly calculable in this 



model fll7 . 



IV. LARGE N DOMAIN WALLS 



As we already mentioned, it is difficult to explicitly construct domain wall solutions in 
the context of SQCD which is a strongly coupled theory. We can circumvent this difficulty 
by considering a theory which is not SQCD itself but is closely related to it, in particular, 
it shares the qualitative features of SQCD such as chiral symmetry breaking and presence 
of domain walls. More precisely, we would like to approximate SQCD by another theory 
which in an appropriate limit reproduces SQCD. We can expect that such an approximating 
theory cannot reproduce (by taking a simple limit of any kind) SQCD at finite A^: the latter 
is a highly complicated model which we cannot hope to be able to solve in a simple manner. 
However, as we discuss in this section, we can write down an approximating model which 
in a certain limit reproduces the large A" SQCD. Using this trick, we, in particular, will be 
able to find the exact solutions for the shapes of the domain walls in the large A" SQCD. 



A. An Approximating Model 

As in section JTj], let us consider Af = 1 gauge theory with SU (N) gauge group and A" 
flavors Qi, Qi (i, j = 1, . . . , N). The non-perturbative superpotential in this theory is given 
by (^|). Next, let us add a tree-level superpotential of the following form: 

(IN 

W tree = XTr(M) - ffjX N+1 ■ (18) 

Here X is a chiral superfield neutral under the SU(N) gauge group. Note that the first 
term in this superpotential corresponds to giving mass |X| to all A" flavors. For \X\ 3> A, at 
energy scales below \X\ we can integrate the heavy flavors out. This leaves us with SQCD 
without matter whose gaugino condensate is controlled by the vev of X. Thus, the effective 
superpotential reads 

W = NX A 2 - -^—X N+1 . (19) 
N + l v ' 

The gaugino condensate is given by 

(AA) = XA 2 . (20) 
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In other words, if \X\ is large compared with A everywhere, then the gaugino condensate in 
this model is a linear function of X. By computing BPS configurations in this model, we 
would therefore be effectively computing the spatial dependence of the gaugino condensate. 

How is this model related to our starting point, namely, SQCD without matter? First, 
let us discuss the vacuum structure of this theory. The F-flatness condition Wx = has iV 
solutions^ (for definiteness here we choose C to be positive): 

(\\) k = x k X 2 , n = (A 2 - N /cy/ N , (21) 

where X k = r]Aexp(2nik/N). Thus, this model has chiral symmetry breaking, namely, 
the Zat subgroup of the .R-parity group is spontaneously broken. There are domain walls 
between different vacua. These domain walls are BPS saturated. This follows from our 
discussion in the previous section, and the fact that this model is nothing but the An model 
when viewed from the two dimensional Af = 2 quantum field theory viewpoint. Thus, this 
approximating model resembles SQCD. However, there are also differences. First, the theory 
contains heavy quarks with masses of order 77 A. In order to be able to approach SQCD, 
we therefore must take the limit 77 —>■ 00. Naively, this might appear to be enough to map 
this model to SQCD. This is, however, not the case. To see this, consider the value of the 
superpotential for a given vacuum: 

N 2 ~ N 2 
W k = X k A 2 = (AA) fc . (22) 

Thus, the relation which would be correct in SQCD, namely, 

W k = N(XX) k , (23) 

is reproduced only in the limit N —>■ 00. (This is required to get the same central charges as 
in SQCD.) So we must also take 't Hooft's large A" limit. In fact, we will see in the following 
that the two limits, namely, 77 — > 00 and N — > 00, are not independent. The precise relation 
between these limits (that is, which combination of 77 and A" needs to be kept fixed) will 
become clear in subsection D. Here we note that in this limit the mass of the field X in 
any given vacuum X = X k is infinitely large, so that this model indeed reproduces large 
A^ SQCD. However, for the domain wall solution interpolating between two distinct vacua 
X = X k and X = X k > there is a zero mode of X localized on the wall. In fact, this zero mode 
appears to be needed for interpreting the corresponding large N domain walls as D-branes 
for the SQCD string. 

B. "Elementary" Domain Walls at Large N 

Next, consider the BPS equation (|TJ) for the above approximating model. (Here and 
in the following we drop the prime and use z instead of z' .) In particular, let us study 



8 There is also the moduli space at X = 0. However, this will have no impact on the following 
discussions since all the configurations of interest here involve I> A. 
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it for the domain wall where at z — > — oo we have X = Xk, and at z — > +00 we have 
X = We have been referring to this wall as an "elementary" wall. First, note that 

exp(— ijk+itk) = i exp(2irik/N + iri/N). Let X = X^Y. Then we can rewrite (|T2| ) in terms 
ofY: 

d z Y* = -texpi-m/N)— [l - Y N ] . (24) 

7] L 

The boundary conditions on Y read: 

Y(z -> -00) = 1 , y(z -> +00) = exp(27ri/iV) . (25) 

Thus, in the large N limit arg(Y) changes by a small amount, namely, from to 2tt/N. It 
is convenient to parametrize Y as follows: 

Y = (1- p/N)exp(i((j) + 7r)/N) . (26) 

Here p is real. The boundary conditions then read: 

p(z -> ±00) = , 0(2 -> ±00) = ±71 . (27) 

Then ( JZty becomes a system of the following first order differential equations: 

9 2 = r[l + exp(-p)cos(0)] , (28) 
3 z p = -Texp(-p) sin(0) . (29) 

Here we are taking the large N limit and only keeping the leading terms. In particular, we 
have taken into account that (1 — p/N) N — > exp(— p) as iV — > 00. Also, we have introduced 

T = N 2 A/ V . (30) 

The above system of differential equations can be integrated. Note that these equations 
do not explicitly contain z. This implies that if <p(z) and p(z) give a solution with the 
appropriate boundary conditions, so will be <fi(z + zq) and p(z + zq) for any constant ^o- 
(This is simply the statement that the system possesses translational invariance in the z 
direction.) There is a solution that has a symmetry with respect to the reflection z — > —z. 
For definiteness, let us focus on this solution. It is given by 

cos(0) = (p - 1) exp(p) , (31) 

Hp) = r\z\ . (32) 



Here 



fPO r -I— i 

^•(p) = / ^ exp(-20 - (1 - 2 2 , (33) 



and po = p(z = 0)(~ 1.278) is the solution of the following equation: 

(po-l)exp(p ) = 1 . (34) 

For illustrative purposes we have plotted the shapes of p(z) and <p(z) corresponding to the 
above exact solution in Fig.l and Fig. 2, respectively. Here we note that the "width" of the 
"elementary" walls is of order 



11 



C. Other Domain Walls at Large iV 



Next, let us analyze the domain walls with k! — k = j, where j is a fixed integer, 
that is, \j\/N — > as iV — > oo. More precisely, we will first consider the cases where 
\j\ ~ 1. In these cases it is convenient to parametrize X as follows. Let X = X^Y where 
Y — (1 — p/N) exp + 7rj)/N). Then we have the following boundary conditions (p is 
real): p{z — > ±oo) = 0, (f)(z — > ±oo) = In the large N limit the BPS equation (|l4l) 

becomes a system of the following first order differential equations: 

d z (j) = T [l - (-l) i exp(-p) cos(0)] , (35) 
d z p=(-iyTexp(-p)sm(<j ) ) . (36) 

It is not difficult to show that for \j\ ^ 1 the only non-trivial solution to this system of 
equations is the one corresponding to \j\ "elementary" domain walls separated by large 
distances (that is, distances much greater than the width of the "elementary" walls). 

Let us understand this fact a bit better. In particular, at first it might appear strange 
that for \j | ^ 1 we do not have any non-trivial solutions other than \j | "elementary" domain 
walls placed at large distances from each other|. Thus, from the two dimensional quantum 
field theory viewpoint we expect to find non-trivial solutions for all values of j. However, 
at large N there are certain simplifications. In particular, consider the energy Ek t k+j of 
the corresponding solitons in the two dimensional quantum field theory. For \j\/N — > 
as iV —>■ oo we have Ek^k+j = 27r|j|/iV = \j\Ek t k+i- This implies that at large N the 
corresponding solitons fall apart into \j\ "elementary" solitons^- In the four dimensional 
language, at large N the domain walls with k' = k+j, \ j\ ^ 1, fall apart into \ j\ "elementary" 
domain walls. In the following we argue that this is consistent with the interpretation of the 
"elementary" domain walls as D-branes for the SQCD string. On the other hand, for finite 
N we have Ek,k+j < \j\Ek,k+i- This is, however, a subleading effect in the large iV limit. We 
therefore conclude that the domain walls with \j\ ^ 1 are stable (in the sense that they do 
not fall apart into "elementary" domain walls) only at finite N. 

To address the above issue in more detail, let us consider the domain walls with k' — k = j 
more carefully. Let us parametrize X as X = XkX . Then the boundary conditions on Y 
read: Y(z — > — oo) = 1, Y(z — > +oo) = exp(27uj / TV) . The BPS equation flUD then reads: 

d z Y* = -i exp(-mj/N) A [l - Y N ] . (37) 

Here we have introduced 

A = NA/r] . (38) 



9 Such a configuration strictly speaking is not a solution. More precisely, for \j\ ^ 1 there are no 
solutions. However, \j\ "elementary" domain walls separated by distances much greater than 
are exponentially close to being a solution of the BPS equation. 

10 More precisely, for finite N the solutions with \j\ ^ 1 can be viewed as "bound states" of 
"elementary" solitons with \j\ = 1 which become infinitely separated as N — > oo. 
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Note that T = NA. 

Suppose there is a BPS solution such that < \Y\ < 1 for z taking values in some finite 
interval. By this me mean that Y N — > as N — > oo for these values of z. Thus, we can try 
the following ansatz (whose self-consistency is to be verified) for solving the BPS equation 



(0): 

3 Z Y* = -iexp(-mj/N)A , z_ < z < z + , (39) 

where it is assumed that Y(z) = 1 for z < Z-, and Y(z) = exp(2mj / N) for z > z + . The 
solution to the above ansatz is given by 

Y(z) = 1 , z < z_ , (40) 
Y(z) = 1 + iAexp(iirj/N)(z - zJ) , z_ < z < z+ , (41) 
Y(z) = exp{2mj/N) , z > z + , (42) 

where z + = z_ + 2L, and 

L = A -1 sin (irj/N) . (43) 

Here z_ is arbitrary. 

Note that there are constraints for the validity of the above solution. First, we must 
consider the "edge" effects. Namely, for z = Z- + e, and z = z + — e, where e C i, we 
have 1 — |y(z)| ~ (e/L) sin 2 (7rj/A^). Now suppose that as iV — > oo we have \j\ ~ N a , 
where < a < 1. Then the assumption Y N —* as iV — > oo breaks down at distances 
e ~ (L/N)N 2{1 ~ a ^ > (L/N) = r -1 . Now recall that T" 1 is the width of the "elementary" 
domain walls. Thus, for \j\ ~ N a , < a < 1, the above solutions would break down due 
to the "edge" effects at distances (from the "edge" ) which are much greater than the width 
of the "elementary" domain walls. This implies that the correct description of the domain 
walls with \j\ ~ N a , < a < 1, appears to be in terms of \j\ "elementary" domain walls 
(separated by large distances) rather than as solitons of the large iV SQCD string. 

Next, consider the domain walls with \j\ ~ N in the large N limit. For such domain 
walls the "edge" effects become important at distances e ~ L/N = T -1 . Thus, the above 
solution for such domain walls can be trusted for all z except for the small regions near the 
"edges". The size of these regions is of order T -1 , that is, the width of the "elementary" 
domain wall. This implies that the domain walls with \j\ ~ iV can be viewed as solitons 
of the SQCD string provided that L is finite and, therefore, the width of the "elementary" 
domain walls (i.e., the width of the D-branes) is vanishing in the large N limit. 

Before we end this subsection, we point out the second restriction on the validity of the 
above solution for the \ j\ ~ iV domain walls interpreted as solitons of the SQCD string. The 
following restriction has to do with the fact the approximating model we have been using is 
strictly speaking valid only if \X\ > A. In the above solution we have |X| = r/A|y|. On the 
other hand, it is not difficult to see that 



min(|y(z)|) = |y(z_ + L)| 



cos 



7TJ 

N 



(44) 



This expression vanishes for j = N/2 (which is possible for even N), and can be arbitrarily 
small for j = (N + j')/2 in the large limit. It is not difficult to see that under these 
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circumstances |AT(z)| < A for z = z_ + L + e, where |e| < L/r]. In the next subsection we 
will see that in the appropriate limit L/r] T~ l . This implies that the above "soliton" 
solutions for \j | ~ A" can be trusted everywhere in the large A" limit up to the "edge" effects 
due to the D-branes, and (possible) "core" effects due to the quarks Qi, Qi becoming too 
light. (These effects, however, are negligible as they affect the solution only in inffnitesimally 
small regions.) 

D. Summary 

The above discussion leads us to the following picture. The width of the soliton solutions 
with \j\ ~ A" is of order L ~ A -1 . In SQCD we expect this to be of order A -1 , where A 
is the dynamically generated scale of SQCD. Note that this is also consistent with the fact 
that L ~ L s , where L s ~ A -1 is the SQCD string length. (The SQCD string tension is 
expected to be T s ~ A 2 .) Then the width of the "elementary" domain walls, interpreted as 
D-branes for the (open) SQCD strings, is of order T _1 = A~ 1 /N ~ L/N, which is vanishing 
for finite L. This is consistent with the (at least naive) expectation that in a weakly coupled 
string theory^ D-branes have vanishing width. 

Another important point is that the domain walls with \j\ ^ 1 appear to fall apart into 
\j\ D-branes (separated by large distances) unless \j\/N is a finite number as A" — > oo. This 
is important for two reasons. Thus, had we found additional domain walls with \j\ ~ N a , 
< a < 1, we would definitely have a puzzle since the tension for these domain walls would 
go as l/Xl +a . Second, it appears that we cannot stack D-branes on top of each other. This 
is consistent with the rest of the above picture, since had we been able to do so, we would 
have to have a moduli space associated with the positions of D-branes. These would have 
to be described by adjoint scalars in the effective field theory. However, we do not expect 
any such adjoint scalars in this model due to M = 1 supersymmetry in 2 + 1 dimensions (in 
the world- volumes of the D-branes - see the next section). So the above conclusions seem 
to be consistent with this fact as well. 

We are now ready to determine the precise limit that we ought to take to recover (the 
large N) SQCD consistent with the above discussions. First, we conclude that we must 
have A ~ A. On the other hand, the gaugino condensate |(AA)fc| = rjA 3 . According to our 
discussion in section ||, we must have |(AA)fc| ~ A^A 3 . Thus, we have two relations: 

NA/r}~A, (45) 
T]A 3 ~ A^A 3 . (46) 

This implies that we must take the following limit: 

T], N -> oo , V /N = C 3 = fixed , A/A ~ C = fixed . (47) 

Here ( is a finite number. (Note that A is fixed to a finite value.) However, we must 
actually have £ ^> 1 (that is, A 3> A) so that the procedure of integrating out the heavy 



The SQCD string theory does not appear to be a critical string theory. The following statement 
is correct for a critical string theory but is a priori far from being obvious for non-critical strings. 
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quarks Qi, Qi is valid. We note that in this limit the above model is a valid description as 
far as the "elementary" domain walls are concerned, that is, none of the assumptions (such 
as, say, validity of integrating out the heavy quarks) on which it was built break down. In 
particular, note that the change in the absolute value of X through the walls corresponding 
to the D-branes is subleading in the large N limit. Thus, the inequality \X\ ^> A always 
holds. Also, the masses of the heavy quarks Qi, Qi (which we have integrated out) are of 
order m ~ \X\ ~ 77A. On the other hand, the width of the D-branes is of order and 
we have m/Y ~ (N/r]) 2 = <C 1. This implies that integrating out these heavy quarks is 
completely legitimate. 

V. ZERO MODES 

The above picture appears to be consistent with the arguments of || as well as various 
discussions in the previous sections. However, if we are to interpret some of the large N 
domain walls as D-branes, we expect (at least naively) that there should be a U(l) gauge 
field living in the world- volume of such a D-brane ||. More precisely, we expect a three 
dimensional M — 1 U(l) gauge supermultiplet localized on the D-brane. Can we identify 
such a supermultiplet in the above picture? 

This is where the zero modes of the chiral superfield X in the above approximating model 
find a natural interpretation. Note that in the above limit the mass of the field X in any 
given vacuum X = Xk is infinitely large. This is consistent with the interpretation of the 
approximating model as reducing to the large N SQCD (where in any given vacuum we 
do not expect additional singlet fields). However, for the solution describing a domain wall 
separating two distinct vacua X = Xk and X = there always exist zero modes of X 
localized on the domain wall corresponding to the Goldstone modes of broken translational 
invariance and supersymmetry. More precisely, since X in the four dimensional M = 1 
supersymmetric theory is a chiral superfield containing a (two-component) chiral fermion 
and a complex boson, the corresponding zero modes (corresponding to a BPS solution that 
breaks half of the supersymmetries) consist of a one-component fermion and a real scalar 
in three dimensions. These form an M = 1 "chiral" supermultiplet in three dimensions. 
However, in three dimensions this is dual to a U(l) vector supermultiplet consisting of a 
vector boson (with one physical degree of freedom) plus its M = 1 superpartner. This is 
due to the fact that a U(l) vector boson in three dimensions can be dualized into a real 
scalar. Thus, for the domain walls corresponding to the D-branes the zero modes of X are 
interpreted (in the dual "magnetic" picture) as the U(l) gauge supermultiplet living in the 
(2 + 1 dimensional) world-volume of the D-brane. 

Let us discuss the zero modes localized on a D-brane in a bit more detail. For simplicity, 
let us confine our attention to the bosonic zero mode. (The fermionic zero mode can be 
obtained by the corresponding supersymmetry transformation.) Thus, consider the domain 
wall separating two vacua X = Xj. and X = X^+i- Let X(z) be the corresponding BPS 
solution. Then the zero mode is given by: 

^(z,x ll ) = d s X(z)<f>(x ll ) . (48) 

Here /1 = 0,1,2. 
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It is not difficult to see that for this D-brane solution the zero mode is localized in the 
region of the size T _1 which is the same as the width of the D-brane. (Recall that in the limit 
we are taking here this width is vanishingly small.) Then according to the above discussion 
the bosonic and fermionic zero modes correspond to the three dimensional M = 1 £7(1) 
gauge supermultiplet localized in the world-volume of the D-brane. 

Here the following remarks are in order. First, note that the above zero modes correspond 
to the translations in the direction transverse to the D-brane. On the other hand, the change 
in the absolute value of X through the corresponding domain wall is subleading in the large N 
limit. Thus, the zero mode is related to the change in the phase of X. The latter is periodic, 
so that the dualization procedure (to a £7(1) gauge multiplet) leads to quantization of the 
corresponding £7(1) charges- 
Second, at low energies we expect a free £7(1) gauge multiplet localized on the D-brane. 
This requires that the dual scalar also be free. This is guaranteed by the fact that here we 
have a Goldstone mode (of broken translational invariance in the z-direction) with purely 
derivative interactions suppressed by the wall tension. (The corresponding interactions 
involving the fermionic component follow by the M = 1 supersymmetry transformation in 
three dimensions.) In particular, at low energies the theory is indeed free. 

Thus, consider the low energy effective field theory corresponding to the zero modes of X 
localized on the D-brane. In particular, let us discuss the bosonic part of the corresponding 
low energy Lagrangian, namely, let us only keep the terms arising through the induced 
metric G^ v on the wall: 

C ~ y/G = det 1 / 2 ( VlMU + d^d^) , (49) 

where G = det(G /Ltjy ). (Here r]^ u is the three dimensional Minkowskian metric, and we 
have normalized the field <fi so that d^cft is dimensionless.) On the other hand, we have the 
following vector-scalar duality transformation in 2 + 1 dimensions (fi, v, X = 0,1,2): 

F^v = e M „A<9> . (50) 

Here F^ u = d^A u — d u A^ is the £7(1) field strength, is the vector potential, and <fi is the 
dual scalar. The Lagrangian in fl4"9|) can be rewritten as 

C ~ det 1 / 2 { Vta/ + F, u ) . (51) 

This is nothing but the Born-Infeld Lagrangian for a £7(1) gauge field (as expected for the 
low energy effective Lagrangian of a D-brane). At low energies it reduces to the free Maxwell 
theory in 2 + 1 dimensions: 

C = -t^V F • (52) 

Here = A 2 F^ llJ has conventional dimension of mass squared (whereas F^ u is dimension- 
less). The £7(1) gauge coupling is given by: 

g D 2 = To/ A 4 , (53) 

where To is the D-brane tension. 



We would like to thank Juan Maldacena for a discussion on this point. 
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VI. DISCUSSION 



Let us briefly summarize the main points in the previous discussions. We consider an 
approximating model which in a particular limit reduces to the large N SQCD. This model 
is constructed by starting from Af = 1 SU (N) gauge theory with N flavors plus a gauge 
singlet. The tree-level superpotential corresponds to the chiral Zn symmetry breaking with 
all the flavors acquiring large masses. There are N inequivalent vacua in this model. In 
a particular limit where the vev of the gauge singlet is large, and N — > oo, this theory 
reduces to the large N SQCD without matter in which we have gaugino condensate. The 
inequivalent vacua then correspond to different phases exp(2iiik/N), k = 0, . . . , N — 1 of 
the gaugino condensate. 

Using our approximating model, we then are able to solve for the exact shapes of the 
domain walls (separating the vacua labeled by k and k! = k + j, j = 1, . . . , N — 1) in the 
large iV SQCD. These exact solutions are obtained in a certain coordinate system which 
depends on the Kahler metric (which we assume to be non-singular). In particular, we find 
that there are two basic types of solutions. The first type corresponds to the "elementary" 
domain walls with \j\ = 1. These are interpreted as D-branes of the corresponding large N 
SQCD string. The tension of these domain walls goes as iV ~ 1/A S (with the appropriate 
choice of the subtraction scheme dictated by the requirement that all domain walls have 
proper interpretation in the SQCD string context, which, in particular, implies that the 
gaugino condensate |(AA)fc| ~ iVA 3 ). The second type of domain walls occur for \j\ ~ N. 
These are interpreted as solitons of the corresponding large N SQCD string. The tension of 
these domain walls goes as iV 2 ~ 1/A 2 . For other values of \j | the corresponding domain walls 
(in the large N limit) appear to be described as configurations composed of \j\ "elementary" 
domain walls (D-branes) separated by distances large compared with the D-brane width. 
(The latter is of order A^/iV — > 0. In contrast, the solitons have finite width of order A -1 .) 

We have also argued that the zero mode (corresponding to the singlet X) localized in 
the world-volume of a D-brane is dual to the M = 1 U(l) vector supermultiplet in 2 + 1 
dimensions. Open SQCD strings should be able to end on these D-branes. Then, these 
D-branes support excitations that behave as if they transform in the fundamental of SU (N) 
||. In other words, as far as the quantum numbers are concerned, the region where the 
SQCD string (i.e., the chromo-electric flux tube) enters the D-brane looks like a state in the 
fundamental of SU(N). This simply follows from the flux conservation. 

The above discussions appear to imply that whatever the SQCD string theory is, it 
must contain D-branes and open strings as well as (expected) closed strings. Here one can 
draw an analogy with Type II string theory which perturbatively contains only the closed 
string states, but must be augmented with D-branes and open string sectors to describe non- 
perturbative phenomena. In fact, in the SQCD string theory it seems to be necessary to 
include the open strings and D-branes (in the large N limit). These states are not solitons of 
the effective closed string theory, but must be introduced or else phenomena such as chiral 
symmetry breaking and appearance of (some of) the domain walls cannot be described. 
At present it is unclear what this SQCD string theory is. However, we do expect that it 
is a non-critical string theory with varying string tension [/J. It would be interesting to 
understand the role of D-branes in such string theories. This might shed additional light on 
the structure of SQCD. 
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Finally, we would like to make the following comment. The particular approximating 
model we used in this paper for our discussion of domain walls in the large N SQCD is 
not the only possible starting point. Thus, one can equally successfully start from, say, 
M = 1 SU(N) gauge theory with one flavor with a large mass. Upon integrating out this 
heavy flavor we obtain SQCD. Then the domain walls can be studied much in the same 
way as we did in this paper. (In fact, this was the setup used in || to study domain walls 
in the case of SU(2).) The resulting domain walls (for large N) are BPS saturated, and 
their shapes in the coordinate system where the corresponding Kahler metric is eliminated 
are related to those obtained in this paper by an appropriate diffeomorphism accompa- 
nied by a (non-linear) change of variables (which changes the Kahler potential). The final 
answer is, however, the same. One should also be able to start from the effective Veneziano- 
Yankielowicz superpotential ||19|| . In this case the corresponding two dimensional M = 2 
quantum field theory reduces to the CP^ -1 model. (See, e.g., |1| for a discussion of this 
model.) Generically, (unless the "glue-ball" field becomes too small) this description should 
also be adequate. Here we should mention that there is a technical problem with using the 
effective Veneziano-Yankielowicz superpotential, namely, a logarithmic branch cut appear- 
ing in the superpotential. It would be interesting to see whether one can obtain domain wall 
shapes in this approach and map them to those discussed in this paper via the appropriate 
diffeomorphism plus change of variables. 
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